Isospin dependent multifragmentation of relativistic projectiles 
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The N/Z dependence of projectile fragmentation at relativistic energies has been studied with the 
ALADIN forward spectrometer at the GSI Schwerionen Synchrotron (SIS). Stable and radioactive 
Sn and La beams with an incident energy of 600 MeV per nucleon have been used in order to 
explore a wide range of isotopic compositions. For the interpretation of the data, calculations 
with the statistical multifragmentation model for a properly chosen ensemble of excited sources 
were performed. The parameters of the ensemble, representing the variety of excited spectator 
nuclei expected in a participant-spectator scenario, are determined empirically by searching for an 
optimum reproduction of the measured fragment charge distributions and correlations. An overall 
very good agreement is obtained. The possible modification of the liquid-drop parameters of the 
fragment description in the hot freeze-out environment is studied, and a significant reduction of 
the symmetry-term coefficient is found necessary to reproduce the mean neutron-to-proton ratios 
{N)/Z and the isoscaling parameters of Z < 10 fragments. The calculations are, furthermore, used 
to address open questions regarding the modification of the surface-term coefficient at freeze-out, 
the N/Z dependence of the nuclear caloric curve, and the isotopic evolution of the spectator system 
between its formation during the initial cascade stage of the reaction and its subsequent breakup. 



PACS numbers: 25.70.Pq,25.70.Mn,21.65.Ef 

I. INTRODUCTION 

Mtiltifragmentation is a universal phenomenon occur- 
ring when a large amount of energy is deposited in a 
nucleus. It has been observed in nearly all types of high- 
energy nuclear reactions induced by hadrons, photons, 
and heavy ions (for reviews, see Refs. [IHlj)- At low 
excitation energies, the produced nuclear system can be 
treated as a compound nucleus which decays via evap- 
oration of light particles or fission. However, at high exci- 
tation energy, possibly accompanied by compression dur- 
ing the initial dynamical stage of the reaction the 
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system will expand to subsaturation densities, thereby 
becoming unstable, and will break up into many frag- 
ments. 

Multifragmentation has been shown to be a fast pro- 
cess with characteristic times around 100 fm/c or less 
(see, e.g., Refs. f9l-[ll|). Nevertheless, as evident from 
numerous analyses of experimental data, a high degree of 
equilibration can be reached in these reactions. In par- 
ticular, chemical equilibrium among the produced frag- 
ments is very likely while kinetic equilibration may not be 
reached [12]. Statistical models were found very suitable 
for describing the measured fragment yields [H y, [13413 ■ 
Even though equilibration may not occur within the indi- 
vidual system during the primary reaction stage the pro- 
cesses are, apparently, of such complexity that dynamical 
constraints beyond the conservation laws do not restrict 
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the statistical population of the available partition space. 

Taking multifragmentation into account is crucially 
important for a correct description of fragment produc- 
tion in high-energy reactions. Depending on the chosen 
projectile, it is responsible for 10 to 50% of the reaction 
cross section. Multifragmentation, furthermore, opens 
a unique possibility for investigating the phase diagram 
of nuclear matter at temperatures T w 3 — 8 MeV and 
densities around p w 0.1 — 0.3po {po ~ 0.15 fm^'^ de- 
notes the nuclear saturation density) . These equilibrium 
conditions coincide with those of the nuclear liquid-gas 
coexistence region. Very similar conditions are realized 
in stellar matter during the expansion phase of super- 
nova explosions [20l - [2^ . The study of multifragmenta- 
tion thus receives a strong astrophysical motivation from 
its relation to stellar evolution and the formation of nu- 
clei, occurring at densities near or below po and over a 
wide range of isotopic asymmetries. The characteristics 
of multifragmentation as a function of the neutron-to- 
proton ratio N/Z of the disintegrating system are, there- 
fore, of particular interest. 

Previous ALADIN experimental data have provided 
extensive information on char ged fragment production in 
multifragmentation reactions |23l. [2^ . In particular, they 
have demonstrated a "rise and fall" of multifragmenta- 
tion with excitation energy [2^ [2^ , and they have shown 
that the temperature remains nearly constant, around 
T ^ 6 MeV, during this process fT2l . [27l [28j , as predicted 
within the statistical approach [29!]. The observed large 
fluctuations of the fragment multiplicity and of the size of 
the largest fragment in the transition region from a com- 
poundlike decay to multifragmentation have been inter- 
preted as manifestations of the nuclear liquid-gas phase 
transition in small systems (see, e.g., Refs. [Ig, I30l433 | 
and references given therein). 

In this article, the study of multifragmentation is ex- 
tended to isotopic effects. New data for projectile frag- 
mentation at 600 MeV/nucleon are presented together 
with their statistical analysis within the framework of 
the statistical multifragmentation model The data 
were obtained during the ALADIN experiment S254 (ssl - 
Issj in which a maximum range in N/Z was explored by 
using primary stable beams as well as secondary radioac- 
tive beams, up to the presently achievable limit on the 
neutron-poor side. The projectile isotopic compositions 
are expected to be essentially preserved during the initial 
stages of the collision, so that the study can fully benefit 
from the asymmetry of available secondary beams. 

The mass of the produced spectator systems varies 
considerably with the impact parameter and is closely 
linked to the transferred excitation energy. The theoret- 
ical analysis is, therefore, based on calculations for a full 
ensemble of spectator systems intended to reflect these 
properties of spectator fragmentation. The parameters 
describing the ensemble of excited nuclei are chosen to 
obtain a best possible reproduction for a set of experi- 
mental observables. We concentrate on observables de- 
scribing the partition degrees of freedom and on the in- 



formation they carry regarding the fragment properties 
at the chemical freeze-out. It is shown that a very satis- 
factory description of the N/Z dependence of multifrag- 
mentation can be obtained within the statistical frame- 
work, including that of the freeze-out temperatures . 
The reduction of the symmetry term coefficient in the 
low-density freeze-out environment, reported for ^^C- 
induced reactions on ^^^'^^"^Sn isotopes is confirmed 
and found necessary to reproduce the measured mean 
neutron-to-proton ratios {N)/Z and the isoscaling pa- 
rameters of intermediate-mass fragments. 



II. EXPERIMENTAL DETAILS 

The ALADIN experiment S254, conducted in 2003 at 
the GSI Schwerionen Synchrotron (SIS), was designed 
to study isotopic effects in projectile fragmentation at 
relativistic energies. Besides stable ^^^Sn, neutron-poor 
secondary Sn and La beams were used to extend the 
range of isotopic compositions beyond that available with 
stable beams alone. The radioactive secondary beams 
were produced at the fragment separator FRS [13] by 
the fragmentation of primary ^^^Nd projectiles of about 
900 MeV/nucleon in a thick berylhum target. The FRS 
was set to select ^^'^La and, in a second part of the exper- 
iment, also ^"^Sn projectiles which were then delivered to 
the ALADIN experimental setup. All beams had a lab- 
oratory energy of 600 MeV/nucleon and were directed 
onto reaction targets consisting of "^'Sn with areal den- 
sity 500 or 1000 mg/cm^. 

For the runs with secondary beams, to reach the nec- 
essary beam intensity of about 10'^ particles/s with the 
smallest possible mass-to-charge ratio A/Z, it was found 
necessary to accept a distribution of neighboring nuclides 
together with the requested ^^"^La or ^"^Sn isotopes. The 
mean compositions of the nominal ^^^La (^''^Sn) beams 
were (Z) = 56.8 (49.7) and {A/Z) = 2.19 (2.16), re- 
spectively [1^. Model studies confirm that these {A/Z) 
values are also representative for the spectator systems 
emerging after the initial stages of the reaction [Sa, . 

The ALADIN experimental setup has been described 
previously [24|. For the present experiment, it has been 
upgraded by constructing a new set of eight proportional 
readout chambers for the TP-MUSIC IV detector and 
by a redesign of all components of the electronic readout 
chain [l^. Digitization of the linear signals was per- 
formed with the aid of 40-MHz sampling ADC's right af- 
ter the preamplifier stage, permitting digital filtering and 
background suppression. The readout planes on either 
side of the TP-MUSIC IV chamber were divided into four 
proportional-counter sections and three ionization-anode 
sections. Each section extended over the full height of 
the chamber, and proportional and ionization sections 
alternated along the beam direction. With this arrange- 
ment, either type of detector could be used to record the 
fragment trajectories. Their horizontal coordinates, i.e. 
in the bending plane of the magnet, were obtained from 
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FIG. 1: Measured cross section as a function of the frag- 
ment atomic number Z determined with the TP-MUSIC IV 
detector for ^'^'^Sn (top panel) and ^"^Sn projectiles (bottom 
panel). Note that the experimental trigger conditions affected 
the yields for Z > 40. 



the measured drift times while the vertical coordinates 
were determined from the positions along the propor- 
tional counters measured with a combination of charge- 
division and pad-readout techniques. 

The threshold for fragment detection and identifica- 
tion was below atomic number Z = 2, obtained with the 
proportional counters. The threshold of the ionization 
sections was higher but they provided a superior Z reso- 
lution for Z > 9. Overall, the resolution AZ < 0.6 (full 
width at half maximum) of the full detector was nearly in- 
dependent of Z up to the projectile Z and permitted the 
individual identification of detected fragments according 
to their atomic number Z. Examples of recorded inclu- 
sive Z spectra for the cases of stable ^^*Sn projectiles 
and radioactive ^^''Sn (nominal) projectiles are shown 
in Fig. [TJ They are given on an absolute scale but the 
recorded cross sections for very large fragments (about 
Z > 40) are underestimated due to inefficiencies of the 
experimental trigger as discussed further below. The fig- 
ure also shows that all elements are produced in mul- 
tifragmentation reactions with an intensity that varies 
rather smoothly with Z. For ^^"^Sn, the recorded inclu- 
sive cross section for the production of fragments with 
3 < Z < 49 is 5.9 b. Charge pickup to Z = 51 ap- 
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FIG. 2: Inclusive mass distributions for fragments with 3 < 
Z < 6 obtained for ^"^Sn projectiles (left panels) and ^^''Sn 
projectiles (right panels). 



pears with 8 mb in the spectrum, a value fitting well 
into the systematics collected for this process |41|. Very 
similar cross sections were obtained for -'^^^Sn projectiles. 
The most significant difference is the lower intensity near 
the atomic number of the projectile in the case of ^°^Sn, 
following from the reduced probability for neutron evap- 
oration of neutron-poor heavy residue nuclei (see below) . 
Traces of the neighboring nuclides in the nominal ^"^Sn 
beams are visible around Z = 50. Smaller differences are 
observed in the range Z < 20 where the odd-even Z stag- 
gering is slightly more pronounced, on the level of 10%, 
in the ^'^^Sn case. 

The magnetic rigidities of detected fragments were de- 
termined by backtracing their trajectories through the 
magnetic field region of the ALADIN magnet, start- 
ing from the vectors determined with the TP-MUSIC 
IV detector. With an iterative procedure, the rigidity 
was varied until the calculated trajectory coincided with 
the measured coordinates of the projectile in the target 
plane, perpendicular to the beam direction. A thin plas- 
tic detector of 200 /xm in thickness, positioned close to 
the target and viewed by photomultipliers at four sides, 
was used for this purpose. The achieved position resolu- 
tion amounted to ~ 1.5 mm, horizontally and vertically, 
over the active area of 4 x 4 cm^ chosen for the secondary 
beam experiments. 
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The fragments, with velocities of approximately that of 
the projectile, can be assumed to be fully stripped of all 
electrons [i^] • Fragment momenta were determined from 
the measured rigidities and atomic numbers Z and used 
to determine fragment masses by including the timc-of- 
flight information given by the ALADIN Time-of-Flight 
Wall [13]. The obtained time-of- flight resolution varied 
with Z, smoothly decreasing from w 250 ps (standard de- 
viation) for lithium fragments to w 100 ps for fragments 
with Z > 10. The start time was derived from the pro- 
jectile arrival time measured with a second thin plastic 
foil of 75 /im thickness, mounted upstream at an angle 
of 45° with respect to the beam direction and viewed by 
two photomultiplier tubes. The resolution obtained for 
the arrival time was « 30 ps (standard deviation). 
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FIG. 3: Reaction cross section dcr/dZbound for ^°^Sn projec- 
tiles as recorded with the interaction trigger (thick solid line) 
in comparison with the expected value obtained with the sub- 
traction technique (dashed line, see text). The solid thin line 
represents, on the same abscissa, the cross section da /dZspcct 
for the production of spectator nuclei with atomic number 
.Aspect according to the geometric participant-spectator model 
of Ref. [H. 

The obtained mass resolution is about 3% for frag- 
ments with Z < 3 (standard deviation), decreases 
to 1.5% for Z > 6, and was found to be the same 
in the experiments with stable and radioactive beams 
(Fig. [2]). Masses are individually resolved for fragments 
with atomic number Z < 10. For optimum resolu- 
tion, correct timing signals from both the front and the 
back layer of the time-of-flight wall had to be required. 
Background-free mass spectra, without spurious peaks 
due to misidentified charges, were obtained by requir- 
ing that the Z values delivered by the TP-MUSIC IV 
and the two layers of the time-of-flight wall all coincided 
within ±0.5 units of the adopted integer value. Because 
of the dead zones between adjacent scintillator slats and 
the reduced charge resolution of the time-of-flight wall 
for larger Z, this condition was met by one-half to one- 
fourth of the fragments with Z < 10. The cross section 
distributions shown in Fig. [5] are, therefore, normalized 
with respect to the recorded Z yields. The significant 



differences of the isotope yields for neutron-poor and 
neutron- rich systems are evident. They are largest for 
the weaker but still identifiable peaks such as, e.g., ^B 
on the neutron-poor or "'^^Be and ^^C on the neutron-rich 
side of the spectra. Isotopically resolved production cross 
sections are shown and discussed in Sec. |Vl 

The acceptance of the ALADIN forward spectrometer, 
in the geometry of the present experiment, was ±10.2° 
for = Z fragments with beam velocity in the horizontal 
direction, i.e. in the bending plane, and ±4.5° in vertical 
direction. At the beam energy of 600 MeV, this corre- 
sponded to an acceptance increasing with Z from about 
90% for projectile fragments with Z = 3to values exceed- 
ing 95% for Z > 6. These values were obtained by fitting 
the measured momentum distributions in transverse di- 
rections with Gaussian functions and by estimating the 
yields lost in the tails. Very similar acceptances have 
been observed in earlier experiments at the same beam 
energy 
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Because of the selective coverage of the projectile- 
spectator decay, the quantity Abound defined as the sum 
of the atomic numbers Zi of all detected fragments with 
Zi > 2 has been chosen as the principal variable for event 
sorting. 2^bound represents approximately the charge of 
the primary spectator system, apart from emitted hydro- 
gen isotopes, and is monotonically correlated with the 
impact parameter of the reaction [2^. The excitation 
energy per nucleon of the spectator system is inversely 
correlated with Zbound 
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Trigger signals for reactions in the target were derived 
from four plastic-scintillator paddles, each 5 mm thick, 
of 8 X 50 cm^ active area, and viewed by photomultiplier 
tubes at both ends. They were positioned approximately 
50 cm downstream from the target at angles outside the 
acceptance of the spectrometer. The condition that at 
least one of them fired was met with nearly 100% proba- 
bility by events with moderate to large charged-particle 
multiplicities. This is demonstrated in Fig. [3] in which 
the rate of interaction triggers obtained in this way is 
compared to the expected event rate. The latter was 
obtained by subtracting event rates recorded with and 
without a reaction target from each other. 

The most notable deviation from the expected rate is 
seen at large Abound at which the expected yields increase 
rapidly while the cross section recorded with the interac- 
tion trigger remains fairly fiat as a function of Abound- A 
rise of the cross section is expected on the basis of the re- 
action geometry as shown in the figure. The smooth thin 
line represents the cross section predicted by the spheri- 
cal participant-spectator model [43] for -'^^^Sn as a func- 
tion of the spectator charge .^spect, plotted on the same 
Z scale. At small Abound, the geometric prediction falls 
considerably below the measured yields because emitted 
hydrogen isotopes are not counted in Abound, but prob- 
ably also because the sharp-edge spherical geometry is 
less realistic for the more central collisions. 

The cross sections dti/dZbound recorded in this way for 
the three cases of ^°^Sn, ^^^Sn, and ^^'^La projectiles are 
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shown in Fig.|4]on a linear scale. Their integrated values 
of 4.3 b (Zbound < 48), 3.9 b (Zbound < 49), and 4.6 b 
(■Z^bound < 55), respectively, compare well with a calcu- 
lated reaction cross section an — 4.8 b (rp = 1.25 fm) 
for Sn + Sn collisions. It confirms that the major part 
of the cross section is recorded with high probability, as 
shown in Fig. [31 Only for very peripheral and extremely 
central events does a bias exist against events with small 
charged-particle multiplicity. 




FIG. 4: Top panel: Triggered reaction cross section as a func- 
tion of Abound for the three projectiles as indicated. Bottom 
panel: The corresponding distributions as obtained from en- 
semble calculations for 500000 reaction events with standard 
parameters (see text). 



III. STATISTICAL APPROACH TO 
MULTIFRAGMENTATION 

A. Reaction scenario and strategy 

In the scenario assumed for the statistical analysis, 
the multifragmentation process is subdivided into (1) a 
dynamical stage leading to the formation of an equili- 
brated nuclear system, (2) the disassembly of the system 
into individual primary fragments, and (3) the subse- 
quent de-excitation of the hot primary fragments. For 
the description of the nonequilibrium stage of periph- 
eral heavy-ion collisions, a variety of models can be used. 
Realistic calculations of ensembles of excited residual 
nuclei that undergo multifragmentation were first per- 
formed with the intranuclear cascade (INC) model |44|. 
More recent examples of hybrid calculations using INC 
or Boltzmann-Uehling-Uhlenbeck (BUU) models for the 
dynamical stage have been reported in Refs. (45l. l46j. 

The dynamical models agree in that the character of 
the dynamical evolution changes after a few rescatterings 
of incident nucleons, when high energy particles ("par- 
ticipants") leave the system. However, the time needed 
for equilibration and transition to the statistical descrip- 
tion is uncertain and model dependent [13, It is 
estimated to be around 100 fm/c for spectator matter 
but the required parameters, i.e., the excitation energies, 
mass numbers, and charges of the predicted equilibrated 
sources, vary significantly even over longer times. To 
evade these ambiguities and other uncertainties related to 
the coupling with a particular dynamical model, alterna- 
tive strategies have been developed [l3 . ll4ll47l . l49l - l53{ and 
are also followed here. With the results of dynamical sim- 
ulations as qualitative guidelines, the exact parameters 
of the thermalized sources are obtained from the analysis 
of the experimental fragmentation data. Usually, only a 
subgroup of observables is required for the determination 
of the source parameters (see, e.g., Refs. [H, |13|)- The 
success of the statistical description is then judged from 
the quality of predictions for other observables and for 
the correlations between them. 



B. Ensemble parameters 



Overall, the cross sections vary rather smoothly with 
■Abound- The drop of the cross section toward the projec- 
tile charge Zq in the neutron-poor cases, while it rises for 
^^''Sn, has already been noted in the Z spectra (Fig. [T]) 
but is observed more clearly on the linear scale (Fig. . 
The fluctuations at small 2'bound are caused by multi- 
alpha events which enhance the even channels. The 
high event rates for Abound = are most likely due to 
5 electrons causing occasional spurious triggers of non- 
interacting projectiles. 



The correlations between mass and excitation energy 
obtained from an INC calculation 1441 and of the geo- 
metric participant-spectator model [43] are shown in the 
top panel of Fig. [S] together with the experimental re- 
suh for i9^Au -f i^^Au at 600 MeV/nucleon [27]. The 
general correlation of an increasing excitation energy per 
nucleon with a decreasing mass of the produced projec- 
tile spectators is a feature common to all cases but the 
values differ considerably. The increase in surface energy, 
solely considered in the geometric model, is much smaller 
than the transferred energies predicted by the INC or ob- 
tained from the experiment. The latter two results are 
fairly consistent with each other but the excitation en- 



6 



ergies were found to be larger than what is needed for 
satisfactory statistical descriptions [l^, 113, IH, A 
considerable part of the accumulated energy is not dis- 
tributed statistically at freeze-out, as evident from the 
fragment kinetic energies which are much higher than 
their thermal expectations js^ ]. 




ensemble of sources - ■ ■ ■ - 
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FIG. 5: Top panel: Mean reduced mass number A/Aq of 
relativistic projectile residuals as a function of their excita- 
tion energy E^/A according to the fireball model of Gosset 
et al. 'isl, calculated for ^°^Sn projectiles, and Intra-Nuclear 
Cascade calculations for ^^^Au-|-C ^4^] (dashed lines), accord- 
ing to the experimental results for ^®^Au nuclei reported by 
Pochodzalla et al. [13] and including the widths in Ex /A 
(dots), and according to the parametrization used in this work 
for 02 = 0.012 and 0.015 MeV"^ (solid lines). Bottom panel: 
Ensemble of hot thermal sources represented in a scatter plot 
of reduced mass number A/Ao versus excitation energy Ex/ A, 
as used in the SMM calculations. The frequency of the indi- 
vidual sources is proportional to the area of the squares. 

In the present work, we follow the procedure used 
by Botvina et al. [l^ which was rather effective in 
describing the multifragmentation of relativistic ^^^Au 
projectiles, including their correlations and dispersions. 
There, the average masses of the equilibrated sources 
As were parameterized as A^/Aq = 1 — ai{E,^/As) — 
a2{Ex/AsY , where is the excitation energy of the 
sources in MeV, and Aq is the projectile mass. The cor- 
relations obtained with parameters oi = 0.001 MeV~^ 



and 02 = 0.015 MeV~^, used in Ref. [11], as weh as 
02 — 0.012 MeV~^ are shown in Fig. [S] It is interesting 
to note that a nearly identical correlation was obtained 
by Desesquelles et al. 51] with a linear backtracing algo- 
rithm based on the SMM and on the fragmentation ob- 
servables measured for ^^^Au -I- Cu at 600 MeV/nucleon. 
As argued in Ref. |47] , a dependence of this kind can be 
considered as consistent with BUU calculations. 

The ensemble of thermal sources adopted for the 
present analysis, for 02 =0.015 MeV~^, is given in the 
bottom panel of Fig. [5] as a weighted distribution in the 
As/Aq versus E* /Ag plane. The widths are empirically 
adjuste d, g uided by the earlier results for ^^'''Au fragmen- 
tation [ij]. The yield distributions as functions of the 
reduced mass A/Aq and excitation energy per nucleon 
were also generated empirically, in a similar manner as 
described in the same reference. Identical parameters 
were used for the three reaction systems. The charges of 
the sources were determined by assuming that the charge- 
to-mass ratio Z/A remains that of the initial projectiles. 
This is expected for the initial reaction phase, as long 
as nucleon-nucleon dynamics dominates, and supported 
by INC and relativistic BUU calculations for similar col- 
lision systems (36l . FbS] . Further arguments in favor of 
this assumption derived from the experimental neutron- 
to-proton ratios of emitted fragments are discussed in 
Sec. ED 



C. Model description 

As the model for multifragmentation, describing the 
reaction stages (2) and (3) defined above, the statisti- 
cal multifragmentation model (SMM, for a review see 
Ref. 0) is used. The SMM assumes statistical equi- 
librium of the excited nuclear system with mass num- 
ber As, charge Zs, and excitation energy (above 
the ground state) within a low-density freeze-out vol- 
ume. All breakup channels (partitions {p}) composed of 
nucleons and excited fragments are considered and the 
conservation of baryon number, electric charge, and en- 
ergy are taken into account. Besides the breakup chan- 
nels, also the compound-nucleus channels are included 
and the competition between all channels is permitted. 
In this way, the SMM covers the conventional evapora- 
tion and fission processes occurring at low excitation en- 
ergy as well as the transition region between the low- 
and high-energy de-excitation regimes. In the thermody- 
namic limit, as demonstrated in Refs. [55l - l57| . the SMM 
is consistent with the nuclear liquid-gas phase transition 
when the liquid phase is represented by an infinite nu- 
clear cluster. 

In the model, light nuclei with mass number A < A 
and charge Z < 2 are treated as elementary stable parti- 
cles with masses and spins taken from the nuclear tables 
("nuclear gas"). Only translational degrees of freedom 
of these particles contribute to the entropy of the sys- 
tem. Fragments with A > 4 are treated as heated nu- 
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clear liquid drops. Their individual free energies Faz are 
parametrized as a sum of the bulk, surface, Coulomb and 
symmetry energy contributions 



Faz - F^ 
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The standard expressions for these terms are F^^ = 
{—Wo — T'^/eo)A, where T is the temperature, the param- 
eter eo is related to the level density, and Wq = 16 MeV 
is the binding energy of infinite nuclear matter; F^^ — 
BoA2/3((t2 - T^)/{T^ + r2))5/4, where Bq = 18 MeV 
is the surface energy coefficient and Tc — 18 MeV is the 
critical temperature of infinite nuclear matter; — 
cZ'^IA^I^, where c = (3/5)(eVro)(l - {plp^f^^) is the 
Coulomb parameter (obtained in the Wigner-Seitz ap- 
proximation) with the charge unit e and = 1.17 fm; 
Ej"^ = 7(A - 2ZflA, where 7 = 25 MeV is the sym- 
metry energy parameter. These parameters are those of 
the Bethe-Weizsacker formula and correspond to the as- 
sumption of isolated fragments with normal density in 
the freeze-out configuration, an assumption found to be 
quite successful in many applications. 

Within the microcanonical treatment 0, [sl] the sta- 
tistical weight of a partition p is calculated as 



(2) 



where Sp is the corresponding entropy, depending on the 
fragments in this partition as well as on the parameters 
of the srotem. In the grand canonical treatment of the 
SMM 59'], after integrating out translational degrees of 
freedom, the mean multiplicity of nuclear fragments with 
A and Z can be written as 



^Naz) 



Yl 

A3 



A^/h 



{Faz -pA- uZ) 



(3) 



Here gAz is the ground-state degeneracy factor of species 

1 /2 

{A,Z), \t = (27rfi,^/m7vr) is the nucleon thermal 
wavelength, and tojv w 939 MeV/c^ is the average nu- 
cleon mass. Vf is the free volume available for the 
translational motion of fragments. The free energy 
Faz = Faz{T,p) of fragments with A> A and Z > 2, 
parametrized as given above, is a function of the temper- 
ature and density. The chemical potentials /i and v are 
found from the mass and charge constraints: 



Y^{Naz)A^As, Y.<^az)Z = Z, 

{A,Z) {A.Z) 



(4) 



For the freeze-out density, one-third the normal nuclear 
density is assumed. This is a standard value, used previ- 
ously in many successful applications and consistent with 
independent ex per imental determinations of the freeze- 
out density [60I l6lj. 

For generating individual multifragmentation events, 
the methods proposed in Refs. 0, [iJI have been used. 
Individual sources are randomly selected from the en- 
semble of excited sources according to their probability 



of occurrence within the ensemble (Fig. [5]) up to a total of 
typically lO'^ to 10^ sources. Their breakup is calculated 
with the SMM by the Monte-Carlo method described in 
Ref. ^ which has the advantage of a high efficiency. In 
particular, at low excitation energies, when partitions of 
primary hot fragments with small total multiplicities M 
dominate, the probabilities of all partitions with M < 3, 
including the compound nucleus formation, are included. 

As a first step, partitions are randomly selected from 
the ensemble of all possible partitions and their prob- 
abilities are determined. For small fragment multiplic- 
ities, the microcanonical sampling according to Eq. [2] 
can be performed. At higher excitation energy, typically 
above E'^ « 4 MeV/nucleon, the masses and charges of 
primary hot fragments are generated according to the 
grand-canonical distribution (Eq. [3]), however, by re- 
quiring the exact mass number and charge conservations 
in each partition. The temperature of a partition is found 
iteratively by adjusting it so that the total energy is con- 
served. It is, therefore, referred to as the microcanoni- 
cal temperature of the partition. After their formation 
in the freeze-out volume, the fragments propagate inde- 
pendently in their mutual Coulomb fields and undergo 
secondary decays. The deexcitation of the hot primary 
fragments proceeds via evaporation, fission, or via Fermi 
breakup for primary fragments with A < 16 [63 [. 

Because the calculated events contain all particles and 
fragments with their masses, charges, and momenta in- 
cluding those produced in secondary decays, they can be 
treated in the same way as experimental events. Sorting 
into event groups can be performed and group averages 
and distributions or correlations can be determined. Be- 
fore comparisons with experimental data are made, the 
resulting partitions of cold fragments are subjected to the 
ALADIN experimental filter by randomly removing frag- 
ments according to the Z-dependent acceptance given in 
Sec. HI 



D. In- medium modification of fragment properties 

The liquid-drop model for the description of the emerg- 
ing fragments is an essential and very successful ingre- 
dient of the SMM 0. It is based on the concept of 
an idealized chemical freeze-out, defined as the moment 
in time at which the nuclear forces have ceased to act 
between fragments. Their action within fragments can 
then be summarized with a liquid-drop description of 
the fragment properties with standard parameters. Only 
the Coulomb energy is reduced due to the presence of 
other fragments, an effect reasonably well accounted for 
within the Wigner-Seitz approximation [2[. However, in 
a more realistic treatment of the breakup stage, the ex- 
cited primary fragments may have to be considered as not 
only excited but also having modified properties due to 
residual interactions with the environment. Their shapes 
and average densities may deviate from the equilibrium 
properties of isolated nuclei [1^ HJ] . Evidence for such 
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conservation are observed throughout this process. 

Possible in-medium modifications of light clusters with 
A < A, treated as elementary particles in the SMM, 
are not considered here. This seems justified because of 
their small primary multiplicities in the niultifragmen- 
tation channels. For example, in the primary partitions 
of a single source ^^"^Sn at = 5 MeV/nucleon, ^'^H 
and ^'^He nuclei appear with mean multiplicities of 0.14, 
0.07, 0.03, and 0.29, respectively. Together they amount 
to a fraction of less than 1.5% of the total mass of the 
system. The majority of light clusters observed in the 
final partitions stem from secondary decays [6^. Mod- 
ifications of light clusters have been intensively studied 
with microscopic approaches, preferentially however in 
temperature-density domains at which larger fragments 
are of reduced importance (see, e.g., Refs. [t^, UM ^^'^ 
references given therein). 



IV. CHARGE CHARACTERISTICS OF 
FRAGMENT PARTITIONS 



10 20 30 40 50 60 



Ensemble normalization 



FIG. 6: Relative yields of fragments for 10-unit-wide inter- 
vals of Abound as obtained with the SMM for the standard 
symmetry-term coefficient 7 = 25 MeV (solid lines) and with 
reduced coefficients 7 = 14 and 8 MeV (dashed and dotted 
lines, respectively) as a function of the fragment Z for the 
case of ^^*Sn projectiles. The centers of the five Zbound in- 
tervals relative to Zo and the scale factors used for displaying 
the yields are indicated in the figure. 



effects has already been obtained from dynamical simu- 
lations [65] and from dynamical or statistical analyses of 
isotopic effects in a variety of reactions [H, [H, [6^ [g^] ■ 
Within the statistical approach, they can be accounted 
for in the fragment free energies by changing the corre- 
sponding liquid-drop parameters. In the following, we 
search for such effects by testing the sensitivity of frag- 
ment observables to the surface and symmetry energy 
terms in comparison with the present data. 

At the last stage of the multifragmentation process, the 
hot primary fragments undergo deexcitation and propa- 
gate in the mutual Coulomb field. The realistic descrip- 
tion of this stage is essential for obtaining reliable final 
fragment yields. To the extent that this can be achieved, 
the hot breakup stage will become accessible, which is 
the preeminent aim of this work. Possible modifications 
of the fragment properties must be taken into account 
in the first deexcitation steps, as the hot fragments are 
still surrounded by other species, while, at the end of 
the evaporation cascade, the standard properties must 
be restored. In the actual calculations, a linear interpo- 
lation between these two limiting cases has been intro- 
duced in the interval of excitation energies between zero 
and E™^ = 1 MeV/nucleon f68|. Energy and momentum 



For event characterization, the total charge bound in 
fragments with Z > 2, 2^bound, has been introduced 
as in previous reports on ALADIN results (see, e.g., 
Refs. p^, M,, 22., 2M)- Small Abound values correspond to 
high excitation energies of the sources that disintegrate 
predominantly into very light clusters ("fall" of multi- 
fragmentation) . Large values of .Abound correspond to low 
excitation energies, at which the decay changes its char- 
acter from evaporationlike/fissionlike processes to multi- 
fragmentation ("rise" of multifragmentation). 

In Fig. m besides the experimental also the theoretical 
event distributions as a function of Abound are shown for 
the three projectiles under investigation, 
and ^'^'^La. The theoretical results are obtained with 
the ensembles introduced in the previous section (Fig. [5]) 
and with standard assumptions regarding the fragment 
properties, i.e. Bq = 18 MeV and 7 = 25 MeV for the 
surface- and symmetry-term coefficients. As .Abound is 
approximately equal to the charge of the spectator sys- 
tem and only slightly affected by its excitation energy, 
the comparison is mainly a test of how well the ensemble 
parameters are chosen. Modifications of the liquid-drop 
parameters have a very small effect on the Abound spec- 
tra. 

The agreement between theory and experiment is over- 
all very satisfactory and sufficient for performing the fol- 
lowing analyses in individual intervals of .Abound- A larger 
disagreement is observed in the region of .Abound > 40 for 
the Sn projectiles and at Abound > 46 for the ^^"^La pro- 
jectile (Fig. 2]). As discussed above, it may be partly 
of experimental origin. The discrepancies at very small 
•Abound < 10 are mainly occurring in channels dominated 
by multi-a events. Here, the larger experimental yields 
may come from more central collisions that are beyond 
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FIG. 7: Mean multiplicity of intermediate-mass fragments, 
(Afimf), as a function of .Zbound as obtained with the SMM 
for the standard symmetry-term coefHcient y — 25 MeV (solid 
lines) and with reduced coefficients 7 = 14 and 8 MeV (dashed 
and dotted lines, respectively), for reactions of ^^*Sn (top) and 
^^*La (bottom) projectiles. 



the statistical treatment of the spectator multifragmenta- 
tion studied here. For these reasons, the following analy- 
sis mainly focuses on the intermediate range of 2'bound ~ 
10-40 (46 in the case of ^"^^La), constituting the major 
part of the rise and fall of nuclear multifragmentation. 

For the quantitative comparison of theory and experi- 
ment, the SMM ensemble calculations were globally nor- 
malized with respect to the measured .Abound cross sec- 
tions in the interval 10 < Abound < 30 for which the 
agreement is best. The obtained factors are 0.00937 
mb, 0.00804 mb, and 0.0106 mb per theoretical event for 
^•^^Sn, ^^"'Sn, and ^^"^La projectiles, respectively, i.e on 
average about 100 events per mb were calculated with 
the SMM. 



B. Sensitivity to liquid-drop parameters 

As a first step, the sensitivity of the fragment-charge 
distributions and correlations to the intrinsic properties 
of these fragments was investigated with the SMM. En- 
semble calculations were performed for different param- 
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FIG. 8: Relative yields of fragments for 10-unit-wide inter- 
vals of Abound as obtained with the SMM for the standard 
surface-term coefficient Bq = 18 MeV (solid lines) and with 
the modified coefficient Bo = 20 MeV (dotted line) as a func- 
tion of the fragment Z for the case of ^^''Sn projectiles. The 
centers of the five .Abound intervals and the scale factors used 
for displaying the relative yields are indicated in the figure. 



eters of the liquid-drop description of the produced frag- 
ments and analyzed as a function of the obtained Abound- 
To permit a more detailed comparison and to include the 
case of ^^^La on the same scale, the reduced ^bound/^Oj 
normalized with respect to the atomic number Zq of the 
projectile, has been used for some observables and subdi- 
vided into five bins (0-0.2, 0.2-0.4, etc.), labeled by their 
central values (0.1, 0.3, etc.) in the legends of the figures. 

For the case of ^^''^Sn, the relative fragment yields as 
a function of Z sorted into these five bins are shown 
in Fig. ini The charge distribution evolves from a so- 
called " U-shaped" distribution, with domination of evap- 
oration of the compound nucleus and asymmetric bi- 
nary decay during the onset of multifragmentation, to 
a rapidly dropping exponential distribution. This evo- 
lution is a characteristic feature of the present type of 
spectator fragmentation [2^ [2^ [t^- The variation of 
the symmetry-term coefficient in the calculations, down 
to one third of its standard value, has a negligible ef- 
fect on the resulting charge distributions except in the 
bin of largest Abound- There the fragment yields in the 
intermediate-mass range which are already small are fur- 
ther reduced if smaller 7 parameters are chosen. 

The same conclusion can be drawn from the results for 
the multiplicities of intermediate-mass fragments (IMF, 
3 < Z < 20) shown in Fig. [7] as functions of Zbound for 
the neutron-rich ^^*Sn and the neutron-poor ^^*La pro- 
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FIG. 9: Mean multiplicity of intermediate-mass fragments, 
(Mimf ) , as a function of Zhound obtained with the SMM for the 
standard surface-term coefficient Bo ~ 18 MeV (solid lines) 
and with the modified coefficient Bo ~ 20 MeV (dotted lines), 
for reactions of ^^*Sn (top) and ^^*La (bottom) projectiles. 



jectiles. The effect of varying the symmetry-term coefB- 
cient is very smah and visible only at the rise of multifrag- 
mentation on the right half of the figure. The faU part is 
expected to be universal because of the auto-correlation 
of the two quantities Mimf and 2'bound lZ3|- The auto- 
correlations dominate in the cases of Abound = 3 and 5 
which can only be reached with partitions containing ex- 
actly one lithium fragment. 

The surface energy of fragments is important for the 
charge yields (tJ- As demonstrated in Figs. [8] and [9l the 
relatively small variation of the surface-energy coefficient 
Bq from 18 to 20 MeV leads to considerable changes in 
the fragment production. The larger surface energy sup- 
presses multifragmentation, leading to smaller IMF yields 
and multiplicities. And, vice versa, decreasing Bq results 
in an intensive breakup of nuclei already at low exci- 
tation energy and to steeper Z spectra. Because of the 
mentioned autocorrelations, the infiuence of Bq decreases 
with decreasing .Abound, as evident from both figures. 
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FIG. 10: Experimental cross sections da/dZ for the fragment 
production following collisions of ^"^Sn projectiles sorted into 
five intervals of Abound /^o with centers as indicated and width 
0.2 (symbols) in comparison with normalized SMM calcula- 
tions (lines). Different scale factors were used for displaying 
the cross sections as indicated. 



C. Data comparison 

The measured cross sections da/dZ for fragment pro- 
duction in reactions initiated by ^''^Sn, ^^*Sn, and 
^'^'^La projectiles are shown in Figs. [T0HT2l The data 
have been sorted into five bins of the reduced bound 
charge ^bound/-^o of width 0.2 spanning the range up to 
•Abound = Zq. The evolution from U-shaped to exponen- 
tial Z spectra, already discussed in the previous section, 
appears without noticeable differences between the three 
cases. 

The results of the SMM calculations, represented by 
the thin lines in the figures, were obtained with the nor- 
malized ensembles (parameter 02 = 0.015 MeV~^, cf. 
Fig. O and with liquid-drop parameters that provide an 
adequate description of the experimental observables and 
trends. The charge observables are mostly sensitive to 
the surface-term coefficient Bq, as discussed in the pre- 
vious section. In order to maintain a consistent compari- 
son, the evidence for a reduced symmetry term coefficient 
7 in the fragmentation channels, extracted from the anal- 
ysis of isotope characteristics in Sec. |Vl has already been 
taken into account by choosing 7 = 14 MeV for the three 
systems. 

For the surface-term coefficient, different values had to 
be chosen in order to obtain equivalent descriptions for 
the neutron-rich and neutron-poor systems. The good 
reproduction of the charge distributions (Figs. [TUHT^ 
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FIG. 11: As Fig.[lO]but for the case of ^^"^Sn projectiles. 
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FIG. 12: As Fig. [10] but for the case of ^^**La projectiles. 



and of the charge correlations shown in Figs. [13] and [14] 
has been achieved with the values Bq = 17.5 MeV for 
^^^Sn and Bq = 19.0 MeV for the neutron-poor ^°''Sn and 
^^^La. They are determined with an accuracy of about 
0.5 MeV for the chosen ensemble of excited sources (cf. 
Figs. [S] and [H]). It should be mentioned here again that 
the comparisons are made with the measured data and 
that the small inefficiencies of the spectrometer for very 
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FIG. 13: Experimental results (top panel) and results of SMM 
calculations (bottom panel) for the mean maximum charge 
(ZmnLx) within an event versus Abound, both normalized with 
respect to Zq, for the three studied projectiles. 



light fragments are taken into account by filtering the 
calculations. If corresponding corrections are applied to 
the data, as done in Ref. , the main effect is a shift of 
the maxima of the mean fragment multiplicity to about 
10% larger values of Abound than shown in Fig. [TU while 
the increase in the maximum multiplicity itself is nearly 
negligible. 

If modifications of the source parameters are simul- 
taneously permitted, the effect of Bq on fragment mul- 
tiplicities can be compensated by a corresponding varia- 
tion of 02. The larger mean multiplicities (Afimf) and the 
shift of the peak of the multiplicity distribution to larger 
Abound observed by reducing the surface term (Fig. [5]) 
can also be obtained by increasing the excitation en- 
ergies of the sources [73|. Equally satisfactory results 
for ^^^Sn were, e.g., obtained with the parameter pair 
a2 = 0.015 MeV-2 and Bq = 17.5 MeV used here as well 
as with a2 = 0.012 MeV'^ and Bq = 19.0 MeV. The 
required N/Z dependence of Bq is independent of this 
ambiguity, at least as long as the ensemble parameters 
remain the same and independent of the projectile N/Z . 

Evidence for a variation of the surface-term coefficient 
with the isotopic composition of the system has previ- 
ously been derived from a t analysis of fragment-charge 
yields y{Z) measured for the fragmentation of ^^^Xe, 

Au, and 238u at 600 MeV/nucleon [7|. Near the 
onset of multifragmentation (.^bound/-^o > 0-6), the r 
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FIG. 14: Comparison of SMM calculations (lines) with the ex- 
perimental results (symbols) for the mean multiplicity (A/imf) 
of intermediate-mass fragments versus Abound for the three 
studied projectiles. 
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FIG. 15: Parameter r obtained from power-law fits to the 
fragment yields yiZ) in the range 5 < ^ < 15 according to 
y{Z) oc Z""^ as a function of Zbound/-Zo- 



additional term 7^ • £-l'^{A - 2ZY/A^ in the liquid-drop 
parametrization. Using the effective A/Z values of the 
two secondary beams and ABq = 1.5 MeV and averaging 
over the two pairs of systems, ^^''Sn vs. ^°''Sn and ^^^Sn 
vs. ^^^La, yields a coefficient 7s = 45 MeV. This is larger 
than parameters of standard mass formulas. That of My- 
ers and Swiatecki, e.g., contains 28.1 and 33.2 MeV for 
the volume and surface coefficients of the symmetry term, 
respectively [zl] . The accuracy of the deduced N/Z de- 
pendence of Bq is not sufficient to draw firm conclusions 
at this point. However, large surface effects are consistent 
with the results for the symmetry term derived from the 
fragment mass yields discussed in the following section. 



parameters obtained from power-law fits according to 
y{Z) oc Z^'^ in the range 5 < Z < 15 are larger for the 
projectiles with smaller N/Z while the opposite is pre- 
dicted by the SMM for constant Bq. There, the variation 
of the N/Z of the studied systems was accompanied by 
a considerable change in mass. The present data confirm 
that the observed inversion of the hierarchy of t with re- 
spect to the predictions with constant B^ is indeed an 
isotopic effect. The r parameters are larger for the two 
neutron-poor systems at large Zbound/-^o (Fig. [15]). For 
the further conclusion that the N/Z dependence of Bq 
may disappear at higher excitations [74| no additional 
evidence is obtained from the present data. In fact, the 
sensitivity of the Z spectra to Bq is seen to diminish at 
smaller .Abound (Fig.|5]) and the fragment multiplicities as 
a function of .Abound become universal (Fig. [5]). 

A variation of the surface term with N/Z is similar 
to a surface dependence of the symmetry energy, i.e. an 



Concluding the analysis of charge observables, it may 
be stated that an excellent overall description of the mea- 
sured yields and correlations is obtained within the sta- 
tistical approach. With the chosen ensemble of excited 
spectator nuclei and after a global normalization, the 
fragment yields and spectra shapes are reproduced on 
an absolute level (Figs. [T0IIT2)) and the correlations agree 
well with the data (Figs. [T^] and [Til) . Deviations oc- 
cur for very large Abound for which the yields of heavy 
fragments are overpredicted and those of light fragments 
{Z < 10) are underpredicted. They should not be over- 
rated, however, because the incomplete trigger efficiency 
in the experiment (Sec. [ll| and the empirical character 
of the yield function of the theoretical source ensembles 
may both be causing them. In fact, the finite values for 
the fragment multiplicity for 2^bound approaching 50 in 
the tin cases seem to indicate a trigger bias favoring the 
more violent collisions. 
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V. ISOTOPE CHARACTERISTICS OF 
PRODUCED FRAGMENTS 

A. Parameter dependence 

The successful description of charge observables pro- 
vides the basis for the investigation of the isotope distri- 
butions and the information contained therein. As in the 
previous section, the parameter dependence is studied 
first. For this purpose, the mean neutron-to-proton ra- 
tios {N) /Z are chosen. SMM calculations of this quantity 
for the two projectiles ^^^Sn and ^^*La and for different 
surface-term coefficients Bq and different symmetry-term 
coefficients 7 are presented in Figs. [16] and [iTl respec- 
tively. 
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FIG. 16: SMM calculations of the mean neutron-to-proton 
ratio (A^) /Z of fragments with 2 < Z < 10 produced in the 
fragmentation of '^■^'*Sn (top panel) and ^^*La (bottom panel) 
projectiles for the interval 0.4 < Abound /^o < 0.6 as a function 
of the fragment Z, performed with 7 = 25 MeV and with two 
surface energy coefficients Bq — 18 and 20 MeV. 
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FIG. 17: SMM calculations of the mean neutron-to-proton 
ratio (N) /Z as in Fig. [I6l but for Bo = 18 MeV and three 
different symmetry energy coefficients 7 = 8, 14, and 25 MeV. 



type of reaction [77] and shown to result from secondary 
decays when the strongly bound even-even N = Z nuclei 
attract a large fraction of the product yields [zl, [t^ . 

The mean neutron numbers are larger for the frag- 
ments of ""^^"^Sn by about half a mass unit, and the dif- 
ference is larger if 7 is reduced (Fig. This differ- 
ence, however, is much smaller than that of the primary 
fragments, which have nearly the same {N)/Z ratio as 
the primary thermal sources [58|. During the secondary 
deexcitation, the neutron richness of fragments will de- 
crease as neutron emission dominates. In the case of 
small 7, this trend is reduced because the mass distribu- 
tion of accessible daughter nuclei is wider and emissions 
of charged particles may more easily compete (681] . 



B. Isotope distributions 



As evident from the figures, the situation is comple- 
mentary to that encountered in the study of the charge 
observables as noted previously [t^. The surface pa- 
rameter B{) has practically no influence on the isotope 
distributions while the symmetry energy strongly affects 
the mean neutron richness of fragments, most noticeably 
for Z >7 . It increases with decreasing 7 because of the 
correspondingly larger widths of the isotope distributions 
at freeze-out. At Z < 6, the mean neutron-to-proton ra- 
tios exhibit the known odd-even effects, typical for this 



Cross sections measured for resolved isotopes and in- 
tegrated over 0.2 < Zbound/^o < 0.8 are presented in 
Figs. [H and [H for Z = 3-6 and 7-10, respectively. The 
major part of the rise and fall of multifragmentation is 
covered by the selected Abound interval (Fig.[T3]). Because 
of the stringent sorting conditions required for produc- 
ing background-free mass spectra, the experimental cross 
sections for isotope production were obtained by normal- 
izing the mass-integrated fragment yields for a given Z 
with respect to the measured elemental cross section in 
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0.2 < Zbound/-^o < 0.8. The dots represent the experimental 
data. Normalized results from calculations with the standard 
and modified parameters are shown by the dashed and solid 
lines, respectively. 



the same .Abound interval (Sec. |TI|. The differences of 
the isotope distributions bet-ween the neutron-poor and 
neutron-rich cases are small but still visible in the loga- 
rithmic representation. In particular, the fragment yields 
on the neutron-rich side are considerably larger in the 
case of -^^"'Sn. 

The SMM calculations sho-wn in the figures -were per- 
formed -with the ensemble parameter 02 = 0.015 MeV~^, 
-with Bq = 17.5 and 19 MeV for the neutron-rich and 
neutron-poor projectiles, respectively, and for t-wo choices 



FIG. 19: Isotope distributions for ^"^Sn (left panels) and 
^^"^Sn (right panels) as in Fig. [18] but for fragments with 
Z = 7-10. 



of 7. For the standard calculations, the symmetry- 
term coefficient 7 = 25 MeV -was used for the hot frag- 
ments at breakup and throughout the secondary deex- 
citation stage. Modified calculations -were performed 
■with 7 = 14 MeV and -with the interpolation procedure 
for the secondary deexcitation a.t < I MeV/nucleon 
(Sec. lIIIj) . The calculated cross sections are globally nor- 
malized for each projectile case as described in Sec. IIVI 
A. 

The comparison with the data sho-ws that the yields of 
neutron-rich isotopes are much better described if the re- 
duced symmetry-energy coefhcient is used. In particular. 
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FIG. 20: Measured mean neutron-to-proton ratios {N)/Z for 
fragments with Z < 10 produced in the multifragmentation 
of ^^'*Sn (top panel), ^''^Sn (middle panel), and ^^''La projec- 
tiles (bottom panel) and for three intervals of Zi,ound/Zo with 
centers as indicated and widths 0.2 (stars). In the calcula- 
tions, performed for the range 0.2 < ^bound/^o < 0.8, the 
surface-energy coefficients _Bo = 18 MeV and the symmetry- 
energy coefficient 7 = 14 MeV were used (open circles). 



the heavier fragments with Z > 7 are rather sensitive to 
the variation of 7, as already evident from Fig. [T71 The 
same conclusion can be drawn from the comparison of 
the mean neutron-to-proton ratios {N)/Z for elements 
of Z < 10 with the experimental results (Fig. I^Ul) . Ev- 
idently, they depend very little on Abound, and the cal- 
culations are, therefore, only shown for the full range 
0.2 < Zbound/2^0 < 0.8. The modified calculations, with 
7 = 14 MeV, agree well with the data, both with re- 
gard to the average magnitude as well as to the odd-even 
structures as a function of Z. The latter, however, are 
slightly more pronounced in the calculations than in the 
experiment for Z < 6 while the opposite is the case for 
larger fragments. This small difference may be connected 
to the discontinuity caused by the mass limit A < 16 set 



for using the Fermi breakup in the secondary-decay part 
of the calculations. 

The ambiguity in the choice of the ensemble and 
surface-term parameters 02 and Bq discussed in Sec. lIVI C 
does not affect this result. In Fig. dTJ calculations are 
shown that were performed with several combinations of 
these parameters, all capable of reproducing the charge 
observables. The observed mean {N)/Z is only reached 
with the reduced symmetry-term coefhcient 7 = 14 MeV. 
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FIG. 21: Mean neutron-to-proton ratio {N)/Z of fragments 
with Z < 10 produced in the fragmentation of ^^'^La. projec- 
tiles in the range 0.4 < .^bound/^o < 0.6. The experimen- 
tal results (stars) are compared with SMM calculations using 
four combinations of symmetry-energy, ensemble, and sur- 
face parameters as indicated. Only the calculations with the 
symmetry-energy coefficient 7 = 14 MeV (circles) are close to 
the data. 

Another important parameter is related to the 
secondary-deexcitation stage of the calculations. For the 
restoration of the standard symmetry-term coefficient, a 
linear interpolation between the initially reduced value 
and 7 = 25 MeV of isolated nuclei in the interval of ex- 
citation energies below E}^^ = 1 MeV/nucleon has been 
adopted (end of Sec. IIII|) . The evaporation times cor- 
responding to this energy are of the order of 10~^^ s 
which seems long, and one may argue that the restora- 
tion should start at a higher initial value fi'i"*, possibly 
closer to 3 MeV/nucleon which would correspond to a 
thermal freeze-out at T « 5 MeV [s^]. Ensemble calcu- 
lations performed with this choice for the ^^^Sn and ^^^La 
systems also provide an excellent description of the data 
(Fig. [22]). It is even superior in the range Z > 6 where 
the odd-even staggering of {N)/Z is much better repro- 
duced than with E'^'^ = 1 MeV/nucleon (cf. Fig.[20l). It 
requires, however, an even larger reduction of the symme- 
try term coefficient 7 to values below 10 MeV. Therefore, 
to the extent that iJ™* — 1 MeV/nucleon has to be con- 
sidered as a lower limit for the beginning restoration of 
the properties of isolated nuclei, the reduced 7 = 14 MeV 
represents an upper limit for the symmetry-term coeffi- 
cient needed to reproduce the observed neutron richness 
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of intermediate-mass fragments. 

Summarizing this section it may be stated that the 
isotopically resolved fragment yields provide evidence 
for a reduced symmetry energy in the hot environ- 
ment at freeze-out, in accordance with previous find- 
ings [H, . Similar observations were made in reaction 
studies performed at intermediate and relativistic ener- 
gies M3LJ3.,M,i8M- 
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FIG. 22: Mean neutron-to-proton ratio {N)/Z of fragments 
produced in the fragmentation of ^^^Sn (top panel) and ^^^La 
(bottom panel) projectiles in the range 0.4 < Zbound/^o < 
0.6. The experimental results (stars) are compared with SMM 
calculations using the three indicated values of the symmetry- 
term coefRcient 7 and _Ei,"* — 3 MeV/nucleon as the start 
value of the interpolation interval for the secondary-decay 
stage of the calculations. With this choice, 7 = 8 MeV (tri- 
angles) gives the best agreement with the data. 



VI. DISCUSSION OF FREEZE-OUT 
PROPERTIES 

A. Source properties 

Several reasonable but also essential assumptions have 
been made in the statistical analysis presented so far. In 
particular, for the N/Z ratios of the considered ensembles 
those of the initial projectiles have been used. At rela- 
tivistic energies, the preequilibrium particles are mainly 
produced in individual nucleon-nucleon collisions and, on 
average, the neutron-to-proton ratios should not change. 
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FIG. 23: Mean neutron-to-proton ratio {N}/Z of fragments 
with Z < 10 produced in the fragmentation of ^'^'^Sn projec- 
tiles in the range 0.4 < Zbound/Zo < 0.6. The experimental 
results (stars) are compared with SMM calculations for the 
nominal N/Z = 1.48 of that system (open circles) and for a 
reduced initial N/Z = 1.33 ratio (dots). The standard liquid- 
drop parameters Bo = 18 and 7 = 25 MeV were used. 



This is contrary to the situation in the intermediate- 
energy range (beam energies less than 100 MeV/nucleon) 
where a considerable isospin diffusion has been identified 
(S^ . Studying this in more detail with INC or BUU cal- 
culations shows that, at relativistic energies, a slight de- 
crease of the neutron richness with increasing excitation 
energy has to be expected. It is not strongly dependent 
on the initial N/Z so that the difference between the iso- 
topic compositions of different residue systems does not 
change by more than a few percent during the initial 
reaction stages [H, HI]. This result is essential for the 
conclusions drawn from the isoscaling study of Le Fevre 
et al. ^B] as well as for the isoscaling analysis presented 
further below. The isoscaling parameter a is roughly 
proportional to the difference of the N/Z ratios of the 
studied systems so that a reduction of the latter, before 
the chemical fragment composition is established, would 
naturally explain the observed decrease of a. 

The situation is different for the analysis of the {N) /Z 
ratios of light fragments presented in the previous section, 
as illustrated in Fig. (23] Two calculations were performed 
for ^^''Sn ensembles, one with the original N/Z — 1.48 
of this projectile and one with a considerably smaller 
neutron-to-proton ratio N/Z =1.33 which is halfway be- 
tween those of the ^^"^Sn and of the neutron-poor ^^^La 
and ^°^Sn nuclei. As expected, the {N)/Z ratios of the 
produced fragments are lower for the smaller initial N/Z 
than for the original choice. To raise them to the exper- 
imental level an even larger reduction of the 7 param- 
eter would be necessary. A larger than assumed varia- 
tion of the isotopic composition of the thermal sources 
at breakup will thus not affect the previously reached 
conclusion. 
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FIG. 24: Apparent temperatures TnoLi (left panel) and TbcLI 
(right panel) as a function of .^bound/-^o for the three reaction 
systems. For clarity, two of the three data sets are slightly 
shifted horizontally; statistical errors are displayed where they 
are larger than the symbol size (from Ref. t35(]). 



B. Freeze-out temperature 

Chemical freeze-out temperatures, as obtained from 
double- isotope-yield ratios measured in this experiment, 
have been presented and discussed in Ref. [35|. The two 
observables shown in Fig. [51] as a function of the normal- 
ized .Abound are the frequently used 

rHcLi = 13.3MeV/ln(2.2^^H/Ji^), (5) 

based on the "^He/^He and ^Li/^Li yield ratios where 
Y denotes the isotopically resolved yields (left panel, 
Ref. ), and 



cLi 



11.3 MeV/ln(1.8^^^^), 



(6) 



the latter being deduced from Li and Be fragment yields 
(right panel, Ref. I28ij). The apparent temperatures as 
given by the formulas are displayed, i.e. without correc- 
tions for secondary decays feeding the ground states of 
these nuclei. Secondary decays are expected to lower the 
apparent values by 10-20% with respect to the actual 
breakup teinperature [27l [28j . 

In Ref. [35j , the near invariance of the measured 
breakup temperatures with the isotopic composition of 
the fragmenting system was shown to be contrary to the 
Hartree-Fock predictions for the limiting temperatures 
of excited compound nuclei 83|. Their global behav- 



ior is, on the other hand, in good qualitative agreement 
with the SMM calculations for ^^^Sn and ^^'^La nuclei of 
Ogul and Botvina [84] . The differences obtained for these 
two cases are negligible in the multifragmentation regime 
and reach a maximum AT « 0.4 MeV in the transition 
region where the equilibrium temperature for the more 
proton-rich ^^^La system is slightly lower. As discussed 
in Ref. [SJ], the difference appears because, in the case 




FIG. 25: Mean microcanonical temperatures calculated for 
the three projectile ensembles with the N/Z dependent sur- 
face coefficients Bq — 17.5 and 19 MeV and with a symmetry 
coefficient 7 = 14 MeV. The experimental result is repre- 
sented by the shaded band that extends from the lowest to 
the highest of the six apparent ThcLI and TbcLI temperatures 
shown in Fig. [24l 



of neutron-rich sources, partitions with a compound-like 
fragment and many free neutrons are still frequent in the 
transition regime. Their temperature is higher than that 
of partitions involving a few charged fragments which are 
more probable for neutron-poor systems. 

The ensemble calculations performed in the present 
work permit a more realistic comparison than those for 
sources with fixed mass and charge of Ref. . The re- 
sult, however, remains the same. As shown in Fig.[25l the 
dependence on N/Z of the mean microcanonical breakup 
temperatures, sorted according to 2'bound, is negligible in 
the regime of multifragmentation. In the most periph- 
eral bin with large .Abound, the temperature of the more 
neutron-rich system is slightly larger, a tendency that is 
also apparent in the experimental isotope temperatures 
(Fig. [24]) . The shaded band in the figure represents the 
full range of the apparent ThcLI and TeoLi temperatures 
which, in the regime of multifragmentation, are close to 
the microscopic temperatures. These results not only 
confirm the usefulness of isotope temperatures for de- 
scribing properties of the chemical freeze-out stage but 
also suggest that the transition from compound decay 
to multifragmentation can be explained statistically with 
the opening of the fragmentation phase space [s^ . 

The calculations were performed with the surface and 
symmetry-energy parameters suggested by the analyses 
of the previous sections, in particular also with 7 = 
14 MeV. Variations of this parameter did not influence 
the results, however [s^. On the absolute scale, the 
microscopic temperature of 6.0 MeV in the center of 
the fragmentation regime ((^bound/^o) — 0-5) is only a 
few percent higher than the apparent isotope tempera- 
tures. The present calculations thus suggest less of a side- 
feeding correction than the up to w 20% derived from 
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a variety of models in Refs. [13, H^. This emphasizes 
once more the quantitative uncertainty inherent in esti- 
mating the effects of secondary decays but also supports 
the conclusion that they are relatively small for the two 
temperature observables chosen here. The calculations 
do not reproduce the accelerated rise of the apparent 
temperature TneLi at very small .^bound/^o ^ 0.2 and 
are still about 1 MeV below the corresponding value of 
TboLi- a contribution to the event class with very small 
■Abound from collisions with more central impact param- 
eters can explain this observation (cf. Sec. HIl. Isotope 
temperatures involving yields of carbon and oxygen iso- 
topes have also been investigated and found to remain 
close to T « 6 MeV, practically independent of Abound 
and in good agreement with the results reported for the 
fragmentation of ^^''Au projectiles [28j . 



C. Isoscaling and the symmetry term 
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FIG. 26: Experimental data (stars) and SMM ensemble cal- 
culations (open symbols) of isoscaling coefficients a extracted 
from fragment yield ratios (3 < Z < 10) for ^^''Sn and "^Sn 
projectiles (top panel) and for ^^'*Sn and ^^"^La projectiles 
(bottom panel) as a function of the reduced bound charge 
Zhou-nd/ Zq. Four different symmetry-term coefficients 7 were 
used in the SMM calculations as indicated. 



Isoscaling concerns the production ratios R21 of frag- 
ments with neutron number N and proton number Z in 
reactions with different isospin asymmetry. It is defined 
as their exponential dependence on N and Z according 
to 

R2i{N,Z)^Y2{N,Z)/Yi{N,Z)<:xexp{N-a+Z-(3) (7) 

with parameters a and /?. Here I2 and Yi denote the 
yields from the more neutron-rich and the more neutron- 
poor reaction system, respectively [s^ HI]. Isoscaling 
has the property that a comparison is made for the same 
isotopes produced by different sources, so that the effects 
of yield fluctuations due to nuclear structure effects may 
be reduced. 
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FIG. 27: Theoretical isoscaling coefficients a calculated for 
three intervals in Zbound/^o and for both hot primary frag- 
ments (open symbols, legend in upper panel) and cold frag- 
ments after secondary deexcitation (closed symbols, legend in 
lower panel). 

The comparison of the measured isotope yields from 
the fragmentation of ^^^Sn and of the two neutron-poor 
systems has confirmed that isoscaling is observed [34, 86] . 
The isoscaling parameter a, determined from the yields 
for 3 < < 10 in different ranges of 2^bound/-^o is shown 
in Fig. [551 It is seen to decrease rapidly as the disin- 
tegration of the spectator systems into fragments and 
light particles increases, as reported previously for the 
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fragmentation of target spectators in reactions of ""^^C on 
"2a24gj^ 3Q0 600 MeV/nucleon Nearly iden- 
tical results are obtained for the isotopic and isobaric 
pairs of reactions. 

In the same figure, the results obtained from SMM 
ensemble calculations with different symmetry-term co- 
efficients 7 and determined for 3 < Z < 10 are com- 
pared with the data. Smaller values of 7 lead to smaller 
isoscaling parameters a because the isotope distributions 
become wider and the variation of the yield ratios with 
N becomes correspondingly smaller. The SMM stan- 
dard value 7 = 25 MeV is applicable only in the bin 
of largest Abound- Smaller values must be chosen to re- 
produce the rapidly decreasing parameter a in the frag- 
mentation regime at smaller Abound- 

The proportionality expected for the dependence of a 
on 7 is demonstrated in Fig. [27l It is fulfilled in good 
approximation for the hot fragments at freeze-out (open 
symbols in the figure). Using the calculated values of 
a for 7 = 14 MeV and the microscopic temperature for 
this interval, T = 6.0 MeV, it is even possible to test the 
widely used formula 



(8) 



where Zi,Ai and ^2j^2 are the atomic and mass num- 
bers of the two thermalized systems. The numerical fac- 
tor deduced from the ensemble calculations is 3.8 while 
the nominal coefficient ^ analytically derived in the zero- 
temperature limit, is 4 (SSf. The validity of Eq.|8]at finite 
temperatures has recently also been confirmed by Souza 
et al. in a theoretical study with the SMM [87]. 

This relation was independently deduced from dynam- 
ical and statistical investigations fsol Issl [89l| and used 
in several analyses of experiments [36, 1661 l67l. l90j . The 
SMM ensemble calculations show that, within some ap- 
proximation, this relation remains valid for the cold frag- 
ment yields as long as 7 is close to its conventional value 
of 25 MeV. As 7 is lowered the resulting a for the asymp- 
totic fragments becomes larger and deviates increasingly 
from the value for hot fragments (Fig. [27)) . During the 
late stages of the secondary decays, the available daugh- 
ter states are increasingly concentrated in the valley of 
stability which narrows the initially wide isotope distri- 
butions. This effect has been demonstrated already and 
discussed in Ref. I36t. 



D. The role of the surface 

The symmetry-term coefficient 7 resulting from the 
isoscaling analysis is shown in Fig. [51] (dots) as a function 
of the sorting variable Abound- For large ^bound/'^o > 0-8, 
the experimental values of a for the isotopic and iso- 
baric pairs of reactions differ (Fig. 1^51) but the predic- 
tion for 7 = 25 MeV is close to the mean. The rapid 
drop for smaller Abound is consistently evident from the 
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FIG. 28: Effective symmetry term coefficient 7 as expected 
from the changing fragment-mass distributions using surface 
and volume symmetry-term coefficients from Refs. [g^, [tsI . 
[9^ . loil j (open diamonds, circles, triangles, and squares, re- 
spectively) in comparison with the coefficient 7 for hot frag- 
ments obtained with the SMM from the N/Z dependence of 
Bo (solid triangles) and from the isoscaling analysis for the 
i07,i24gj_^ pair of reactions (dots). 



two data sets (Fig. [261) . For comparison, four predic- 
tions are shown in Fig. I28[ obtained with rather different 
approaches but all containing the effect of the surface- 
symmetry term whose importance increases for the lower- 
mass fragments [9T'| . The coefficients of the mass formula 
of Myers and Swiatecki [75i| are adapted to ground-state 
masses, and values close to them have been used in other 
studies [9^. From the energies of isobaric analog states, 
a relation between the volume and surface capacitances 
of nuclei for absorbin g as ymmetry N — Z was derived by 
Danielewicz and Lee i93|, while Kolomietz and Sanzhur 
have used a variational approach using Skyrme forces to 
derive equilibrium values for the volume symmetry term 
with surface and curvature corrections for nuclei along 
the /3-stability line [sJ]. Theoretical surface and volume 
terms were also obtained by Ono et al. ^65'| from fitting 
the ground-state binding energies calculated with the an- 
tisymmetrized molecular dynamics (AMD) model. 

With these parametrizations, an effective symmetry 
energy averaged over the set of partitions was calculated 
for the five bins in Zbound after the experimental Z distri- 
butions had been converted to mass distributions using 
the projectile N/Z. The obtained results show similar 
trends. The smaller fragments produced at higher excita- 
tions cause the effective mean symmetry term to decrease 
with decreasing Zbound in all four cases but at a slower 
rate than that resulting from the isoscaling analysis of 
the experimental yield ratios. 
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The same procedure was, in addition, also performed 
with the surface-symmetry coefficient 7^ — 45 MeV de- 
rived from the N/Z dependence of the surface-term coef- 
ficient Bq (Sec lIVj) . combined with the volume-symmetry 
coefficient j„ = 28 MeV of the mass formula of Myers and 
Swiatecki [75|. As expected, the drop with decreasing 
■Abound is faster because of the larger surface-symmetry 
coefficient and close to that deduced from the isoscaling 
analysis. Within their larger uncertainty, the conclusions 
drawn from the analysis of the charge observables are 
thus consistent with those obtained from the isoscaling 
analysis of the measured isotope distributions. Increased 
surface effects are expected if deformations and exotic 
shapes are present among the smaller fragments at the 
chemical freeze-out. 

The symmetry-term coefficients of the AMD ground- 
state nuclei follow the trend exhibited by the other mass 
formulas. However, for excited nuclei at the breakup 
stage, a reduced volume term and an isospin-independent 
surface term for the larger IMF's were found in the AMD 
study [65l|. From this result and from the t analysis of 
fragment charge yields [7J] one may conclude that the 
strong decrease of 7 with decreasing .Abound is unlikely 
to result from a surface effect alone. In fact, the surface 
term should decrease also for fragments in close vicinity 
of other nuclei because of residual interactions. The re- 
duced symmetry term will then have to originate from 
a decrease of the overall matter density caused by the 
expansion of the system [69l . losj and leading to expanded 
or deformed shapes of fragments as they are produced. 
These effects should be considered in the AMD model 
but were, perhaps, not tested over the full range of pos- 
sible excitations in Ref. [6^1 • There, the temperature at 
breakup was estimated to be less than 4 MeV which cor- 
responds to the bin of largest Zbound here (Fig. [M]) for 
which the deduced 7 is comparable to the AMD result. 
The stronger reduction is observed at the higher temper- 
atures at smaller Abound (Fig- [M]) . 

On the quantitative level, different degrees of reduction 
are deduced with the different methods. The lowest val- 
ues for 7 are obtained from the isoscaling analysis. The 
7 « 10 MeV obtained for the center of the fragmentation 
regime, i.e. .^bound/-^o ~ 0-5) is to be compared to 13 
MeV, corresponding to ABq ~ 1.5 MeV but carrying a 
large uncertainty, and to 7 = 14 MeV which provided 
the best global agreement with the charge and isotope 
observables studied in Sees. IIVI and |Vl An answer may 
come from the fact that isoscaling is more sensitive to the 
widths of the isotope distributions and thus may give a 
larger weight to their tails. These quantities, on the other 
hand, may more strongly depend on the realization cho- 
sen for the secondary decay stage of the reaction. As 
discussed above, values below 14 MeV are favored by the 
odd-even structure of {N)/Z for Z > 6 (Fig.B^. 

The global 7 = 14 MeV seems close to the effective 
coefficient 7 = 17 MeV obtained from the parametriza- 
tions of Myers and Swiatecki or Danielewicz and Lee, 
including only surface effects of isolated nuclei (Fig. [251) . 
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FIG. 29: Mean neutron-to-proton ratio {N)/Z of fragments 
with Z < 10 produced in the fragmentation of '^^''Sn (top 
panel) and ^^"^La (bottom panel) projectiles in the range 
0.4 < Zbound/^o < 0.6. The experimental results (stars) are 
compared with SMM calculations with the Myers-Swiatecki 
parametrization [t^ for the liquid-drop description of frag- 
ments and for two ensembles as indicated (open symbols). 

However, the difference is significant. SMM calculations 
using the coefficients of Myers and Swiatecki for the frag- 
ment description throughout, i.e., at breakup as well as 
during the secondary decay stages, lead to significantly 
smaller {N)/Z, as shown for the ^^"^Sn and ^^^La. sys- 
tems in Fig. 1291 They are not sufficient to reproduce 
the observed neutron richness of light fragments. When 
searching for possible reasons for these large differences it 
was found that the reduction of 7 lowers the Q- value for 
proton emission while this is not the case with the Myers- 
Swiatecki parameters. The relative probability for proton 
emission is thus enhanced, leading to more neutron-rich 
daughter nuclei. 

VII. CONCLUSIONS 

New results from a systematic study of the N/Z depen- 
dence of projectile fragmentation at an incident energy 
of 600 MeV/nucleon have been presented. By employing 
radioactive secondary beams, a large range of isotopic 
compositions has been explored, up to the present tech- 
nical limit. Global fragmentation observables were found 
to depend weakly on the projectile N/Z while substan- 
tial differences are observed for the isotope distributions 
of light fragments {Z < 10). 
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The significance of the measured charge and mass 
yields and correlations was investigated with the SMM 
by performing ensemble calculations accounting for the 
mass-energy correlation of the produced spectator nuclei. 
The ensemble parameters were determined empirically by 
searching for an optimum reproduction of the measured 
fragment Z yields and correlations. The sensitivity of 
the liquid-drop parameters of the fragment description 
in the hot freeze-out environment was studied and, as 
a main result, a significant reduction of the symmetry- 
term coefficient 7 was found necessary to reproduce the 
mean neutron-to-proton ratios (TV) jZ and the isoscaling 
parameters of Z < 10 fragments. With 7 = 14 MeV 
a globally good reproduction was obtained while even 
lower values are suggested by some observations. Stud- 
ied as a function of the impact parameter, 7 is found to 
decrease rapidly as the multiplicity of fragments and light 
particles from the disintegration of the produced specta- 
tor systems increases. The reduction is larger than that 
expected from the decreasing fragment size and the cor- 
respondingly larger surface effects as modeled with stan- 
dard mass formulas. It is interpreted as resulting from 
the overall reduced density of the system at breakup and 
the modified fragment properties in the hot environment. 
Deformations and reduced fragment densities at chemical 
freeze-out contribute to producing the observed yields of 
neutron-rich fragments which, according to the present 
investigation with the SMM, are particularly sensitive to 
the widths of the isotope distributions and thus to the 
symmetry-energy term at freeze-out. 

The reproduction of the experimental data with the 
SMM ensemble calculations is found to be very satisfac- 
tory. A global normalization on the level of the .^bound- 
sorted event distributions was found sufficient for reach- 
ing an overall quantitative description of the Z and A 
resolved fragment yields. The calculations have, there- 
fore, also been used to address open questions regarding 
the isotopic compositions and temperatures of the sys- 
tem at the freeze-out stage. The N jZ independence of 
the breakup temperatures in the fragmentation regime, 
obtained from double-isotope yields of light fragments, 
is quantitatively reproduced by the microscopic SMM 



temperature. The good agreement suggests that the 
opening phase space is sufficient to initiate the tran- 
sition from predominantly residue formation to multi- 
fragment breakups as the transferred energy increases. 
Changes of the isotopic evolution of the spectator sys- 
tem between its formation during the initial nucleon- 
nucleon cascades and its subsequent breakup have been 
ruled out as alternative explanations for the isotopic frag- 
ment compositions. A reduced difference of the system 
N jZ at breakup will result in smaller isoscaling parame- 
ters but, at the same time, require even larger reductions 
of the symmetry-term coefficient in order to reproduce 
the mean (iV) jZ of intermediate-mass fragments. 

It is finally emphasized that the simultaneous analysis 
of several observables obtained from a nearly complete 
coverage of the spectator kinematic regime has been im- 
portant for reaching the presented conclusions. The de- 
duced modifications of hot fragments in the low-density 
freeze-out environment and in the presence of other frag- 
ments and nuclei can be expected to influence the mod- 
eling of astrophysical processes. 
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